Critical dynamics, duality, and the exact dynamic exponent in extreme type II 

superconductors 
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The critical dynamics of superconductors is studied using renormalization group and duality 
arguments. We show that in extreme type II superconductors the dynamic critical exponent is 
given exactly by z = 3/2. This result does not rely on the widely used models of critical dynamics. 
Instead, it is shown that z = 3/2 follows from the duality between the extreme type II superconductor 
and a model with a critically fluctuating gauge field. Our result is in agreement with Monte Carlo 
simulations and at least one experiment. 
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I. INTRODUCTION 



The high-T c cuprate superconductors have very large 
values of the Ginzburg parameter n. In a material 
like YBa2Cu307_a (YBCO) at optimal doping, we have 
k ~ 100. It is therefore a good approximation to as- 
sume that these materials are extreme type II supercon- 
ductors. By extreme type II superconductor we mean 
k — ► oo . In such a regime it is expected that the static 
critical properties at zero external field are the same as in 
superfiuid 4 Hc,i This expectation is confirmed by exper- 
iments measuring specific heali^ and penetration depth 3 
in bulk samples of YBCO. Thus, as far as static criti- 
cal phenomena are concerned, there is no doubt that the 
critical behavior of bulk YBCO is governed by the three- 
dimensional XY universality class. On the basis of such 
a consensus, we might also expect that the dynamical 
universality class is the same as in superfiuid 4 He. Un- 
fortunately, we are far from reaching a consensus on the 
dynamical universality class of YBCO, or more generally, 
of any high-T c cuprate. What is interesting here is that 
the lack of consensus comes both from the theoretical and 
experimental sides. The theoretical debate tries to estab- 
lish whether the dynamical universality class corresponds 
to model A or model F dynamics. 4 Model A dynamics is 
purely relaxational and gives the value z w 2.015 for the 
dynamic critical exponent in three dimensions. Model F, 
on the other hand, features a conserved density coupled 
to a spin- wave mode and gives the exact value z = d/2 
for d G (2,4]. For the dynamical universality class of 
extreme type II superconductors, Monte Carlo simula- 
tions give z ~ 3/2j££ which would be consistent with 
model F dynamics and therefore with superfiuid 4 He dy- 
namical universality class. However, in a recent letter 
Agi and Goldenfeldi claim that the correct result from 
a lattice model should be instead z « 2, i.e., consistent 
with model A dynamics. From the experimental side it 
is found either large values for z, typically in the range 
z ~ 2.3 — 3.0, 8 « either values consistent with model F 



dynamics^ The arguments of Ref. were further dis- 
cussed in two recent comments^ 

In principle model F dynamics is not compatible with 
superconductor critical dynamics, since screening ef- 
fects tend to suppress the spin- wave modeil^ However, 
model A is also inappropriate to study the critical dy- 
namics in superconductors. Indeed, model A does not 
give a gauge-independent result for z in the magnetic 
critical fluctuation (MCF) regime. A technically cor- 
rect analysis should consider the extreme type II regime 
as a limit of the full MCF regime which has a genuine 
local gauge symmetry. It turns out that in a gauge- 
invariant theory the only operators with a nonzero ex- 
pectation value are the gauge-invariant ones*^ For in- 
stance, the AC conductivity is such a gauge-invariant 
quantity. Therefore, its scaling behavior is necessarily 
gauge-independent. Since evaluation of z through model 
A in the MCF regime gives a gauge-dependent result, we 
conclude that model A is also not the right option. 

In this paper we use renormalization group (RG) and 
duality arguments to determine the dynamical universal- 
ity class of extreme type II superconductors. We will 
establish that in three dimensions z — 3/2 exactly. This 
result will be obtained by combining exact scaling ar- 
guments with exact duality results. While the scaling 
arguments are generally valid in d G (2,4], the duality 
arguments will be only valid at d = 3. The result will be 
obtained through the following strategy. In Section II wc 
use the RG to obtain exact scaling relations for the pen- 
etration depth and for the AC conductivity. This will be 
done in both the XY (extreme type II limit) and MCF 
regime. While much of the steps in this part of deriva- 
tion are known, it is important to review this approach 
here to emphasize the intimate relationship between the 
scaling of the penetration depth and the one of the AC 
conductivity. The important point is that scaling behav- 
ior in the XY regime is different from the one in the 
MCF regime. In Section III we discuss critical dynamics 
from the perspective of the exact duality between an ex- 
treme type II superconductor at zero field and a model 
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exactly equivalent to the Ginzburg-Landau (GL) model 
in the MCF regime. Indeed, the dual model features a 
fluctuating vector potential coupled to a complex field, 
dubbed disorder parameter field as opposed to the order 
parameter field of the original model. The derived re- 
sults for the GL model in the MCF regime will then be 
applied to the dual model and duality relations between 
the currents will be used to establish that z = 3/2. 



II. SCALING AND DYNAMICS IN THE 
GINZBURG-LANDAU MODEL 

In order to make the paper self-contained, we review in 
this Section the basic scaling properties of the GL model 
in both the static and dynamic critical regimes. 



A. Static scaling behavior 

Let us consider the bare Hamiltonian of the GL model, 



ones. From this it can be shown that the following exact 
evolution equation for the renormalized vector potential 
mass holds 



(3) 



where 7a(/, g) — ndhiZ 'a/ 'dp and (3 g = \idgjd[i. Eq. 
© implies that near the phase transition 



(d-2+ VA )/2 



(4) 



where r\A — l\ is the anomalous dimension of the vector 
potential. In the MCF regime r] A = 4 - rf^ 18 . 20 . 24 and 



we obtain that v 



20 



where v' and v are the penetra- 



tion depth and correlation length exponents, respectively. 
This result was confirmed by Monte Carlo simulations. 22 
In the XY regime, on the other hand, the penetration 
depth exponent is given by«ii 



v(d-2) 



(5) 



H = i(VxA ) 2 



KV-zgoAo^op+^I^P + yl^ol 4 , 



(1) 

where qo = 2eo is the charge of the Cooper pair and 
/1q oc r, with t — (T — T c )/T c . In our notation the 
zero subindex denotes bare quantities while in renormal- 
ized quantities the zero subindices are absent. The bare 
Ginzburg parameter is kq — Ao/£o = {u-o/^Qo) 1 ■> where 
Ao and £o = Mo" 1 are the bare penetration depth and cor- 
relation length, respectively. We can rewrite the above 
Hamiltonian in terms of renormalized quantities as fol- 
lows: 



H 



?A(VxAf+z^\(v^ q A)ijA 2 +z^ 2 M 2 +^rM 4 



The renormalized correlation length is given by £ = fi 1 . 
From the Ward identities it follows that the renormalized 
charge squared is given by q 2 = ZaOq- The dimensionless 
couplings are / = fi d ~ 4 q 2 and g = [i d ~ 4 Z 2 uo/Z u . The 
fixed point structure is well known but cannot be com- 
pletely obtained by perturbative means. Fixed points as- 
sociated to nonzero charge, /* ^ 0, are non-perturbative 
but their existence in the flow diagram is now well 
establishedi 15 i 16 . 17 i 18 i 19 i 20 . 22 i 23 i 24 i 25 The infrared stable 
charged fixed point governs the MCF regime while the 
extreme type II or XY regime is governed by the un- 
charged XY fixed points. 

The renormalized Ginzburg parameter is given by k — 
fi/HAi where \i A = A -1 is the renormalized vector poten- 
tial mass generated by the Anderson-Higgs mechanism. 
Due to the Ward identities this can also be written as 
k — (g/2/) 1 / 2 , and therefore the renormalized Ginzburg 
parameter has the same form as the bare one, with the 
bare coupling constants replaced by the renormalized 



where v is the correlation length exponent. At d = 3 
we have v « 2/3 and v' w 1/3. The result v' 1/3 is 
confirmed experimentally in high quality single crystals 
of YBCO£ Note that in the XY regime kq — > oo while 



B. Dynamic scaling behavior 

The AC conductivity is given by 



cr(w) 



q 2 K{-iu) 



(6) 



where K (— iu>) is obtained from the current-current cor- 
relation function at zero momentum^ i.e., K(—ioj) — 
ftm\ p \_> K(-iu),p) where K(-iu,p) = K ^(~ iuJ iP) 
with 



K^i-iu^p) = (\M 2 )S„u - -5<J p (w,p) • Ju(-u,p)), (7) 

and J M (w,p) is the Fourier transform of the supercon- 
ducting current 



(8) 



The function K^ v (~iu, p) is purely transverse^ 7 - The su- 
pcrfluid density p s is given by lim^^o K(—iu>) and thus, 
in virtue of the Josephson relation, we obtain 



lim K(-iuj) ~ (-iio) (d - 2 ^ z . 

f-L—>-0 



(9) 
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Since q 2 ~ fi nA , we obtain from Eq. © the behavior 

<r(w)| r =T B ~ (-^)( d - 2 ^ z +^)/ z . (10) 
Therefore, in the XY universality class we hav o 11 ' 12 



<7(W)| T =T. ~ H^- 2 "^, 

while in the MCF regime we obtain 



a(u))\ T =T a 



(-iui) 



(2-z)/z 



(11) 



(12) 



From Eq. © we obtain that below T c the low fre- 
quency behavior of the AC conductivity is 



(13) 



In Refiii this is written symply as a(oj) w p s /(—ioj), 
i.e., the charge is not shown explicitly since only the XY 
regime was considered and in this case the charge does 
not fluctuate. Taking the low frequency behavior ljl3|l 
into account, we obtain from Eqs. © and |JBJ|, 



/'- 



da{uj) 



d - 2 — z + j a 



ft i , , 
1 cr(w), 



V 5 
which implies that near the phase transition, 



(14) 



d— 2— z+t)A 



<j(u>) ~ /Lt 



Since in the XY universality class t\a 
from Eq. I|15(l the well known scaling 11 



cr(u>) ~ /i' 



d-2- 



(15) 

0, we recover 
(16) 



Note that Fisher e£ ali^ need to assume the Josephson 
relation p s ~ ^ rf_2 to derive the X7 scaling of the AC 
conductivity. Within our approach the more general scal- 
ing relation l|15|l follows from Eq. © and the XY scal- 
ing emerges as a particular case. In the MCF regime we 
obtain 29 ' 30 



<j(oj) ~ /i 2 



(17) 



III. DUALITY AND CRITICAL DYNAMICS 

A. Duality and disorder field theory 

In the extreme type II limit kq — > oo and we have essen- 
tially a superfluid model at zero field, i.e., the correspond- 
ing Hamiltonian is the same as in Eq. (TJ with A = 0. 



The lattice version of this model in the London limit is ex- 
actly dual to the so called "frozen" superconductor 32,33 . 
Starting from the Villain form of the XY model we ob- 
tain, after dualizing it, the following lattice model Hamil- 
tonian: 



H = J2 



i 

2K 



(V x h,) 2 - 2m M t • h, 



(18) 



where K is the bare superfluid stiffness, a; M € (—00,00) 
and Mi^ is an integer link variable satisfying the con- 
straint V • M; = 0. The lattice derivative is defined as 
usual, V M /; = fi+jx — The link variables play the 
role of vortex currents and the zero lattice divergence 
constraint means that only closed vortex loops should be 
taken into account. Integration over a/ gives a long range 
interaction between the link variables. The link variables 
will interact through a potential V(ri — r m ) behaving at 
large distances like V(yi — r m ) ~ l/|r; — r m |. At short 
distances the potential is divergent. This short distance 
divergence can be regularized by adding to the Hamilto- 
nian (|18(l a core energy term (eo/2) J2i Mf. Writing the 
constraint V • M; = using the integral representation 
of the Kronecker delta and performing the sum over M/ 
using the Poisson formula we arrive at the Hamiltonian 



1 

2^ 



(V0j - 2ttN z - 27TV K h 



(19) 



where we have rescaled hj. The dual lattice Hamilto- 
nian l|19|) has exactly the same form as the Hamiltonian 
for a Villain lattice superconductor. Note that 2-k\[K 
plays the role of the charge. The sum over the integers 
Nj can be converted into a disorder field theory (DFT)JA 
which has precisely the same form as the original GL 
model in Eq. QJ, except that the physical properties of 
the fields have changed. The electromagnetic vector po- 
tential Ao is replaced by the gauge field ho describing 
vortices, and the charge go becomes the Biot-Savart-type 
coupling strength between vortices 2tt^/p2, where p° is 
the bare superfluid density. An important aspect of du- 
ality is that the disorder field theory that comes out of 
it has an "inverted" temperature axisji^ i.e., the broken 
symmetry phase of the disorder field theory corresponds 
to the symmetric phase of the original theory and vice- 
versa. This is actually the meaning of the expression 
"disorder field" used in this paper: instead of having an 
order 'parameter as in the original model, the dual model 
has a disorder parameter. Both models describe the same 
physics and have of course the same critical temperature 
T c . Since the Ward identities for either theory imply 
that the critical singularities are the same irrespective of 
whether T c is approached from above or from below, we 
can use the same scale /1 to study the scaling behavior in 
both models. 23 ^Si The critical exponents a and v of the 
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dual model are the same as in the original model. This is 
because the dual model give the same free energy of the 
original model, up to non-singular terms. Therefore, the 
exponent appearing in the scaling of the singular part of 
the free energy, a, is the same in both models. The hy- 
perscaling relation then implies that v is also the same 
in both models. 

The continuum dual model for superconductors was in- 
troduced in Ref. and further discussed in Ref. (see 
also Ref. l36|) . Such a continuum dual theory represents a 
generalization of the London model. In the extreme type 
II limit and zero external magnetic field, it is obtained 
from the continuum limit of Eq. i|19|) , 



7i-i(Vxh o ) 2 + |(V-jgoho)0o| 2 + m2|0 o | 2 + ^|0 o | 4 , 



where the dual bare charge is given by 



(20) 



Thus, by taking the curl of Eq. (|22() . we obtain a dynam- 
ical equation for the vortex current: 



dJ v 2 fiA,o dB 

— -T oPs V J v -— — , 



(24) 



where B is the macroscopic magnetic induction field. In 
Fourier space and at zero momentum we obtain 



3 v (u>) 



fJ-A,o 
' 2ir 



B(w). 



(25) 



Thus, the above simple analysis already shows us that the 
introduction of dynamics in the duality approach leads 
to interesting consequences. In the original model linear 
response theory relates the current to the electric field in 
Fourier space as 



go = =2ttVp^ (21) 

q 

The Hamiltonian 1)20(1 has the same form as the GL 
Hamiltonian in Eq. In Eq. I|20|) the gauge field ho 
is minimally coupled to the already mentioned disorder 
field 4> . We should note that this disorder GL-like theory 
is valid strictly in d = 3. Indeed, if we try to extrapolate 
to the range of dimensionalities as in scaling relations 
of Section II, we would obtain that the renormalized su- 
perfluid density, as the "charge" of the DFT, scales like 
p s ~ jtx 4 , which would agree with Josephson's relation 
only for d = 3. 

B. Critical dynamics of the dual model 

In the duality transformations in the lattice only the 
phase of the order paramerer plays a roleiSL^i Actually, 
often duality arguments are worked out directly in the 
continuum, by relating the gradient of the phase to a 
"magnetic field" that couples to the vortex loops within 
the dual modeli2S4£ In the case of critical dynamics it 
is simpler to follow a similar approach in order to derive 
the main results. 

The superfluid velocity w s satisfies the dynamical 
equation 2 ^ 



J(w,p)=ff(w,p)E(w,p). (26) 

With respect to the dynamics of the system, the sta- 
tistical mechanics duality implies also electric-magnetic 
duality, in which case the electric field is replaced by the 
(true) magnetic field in a linear response theory for the 
vortex current. 

Generally, the current of the original GL model is re- 
lated to the vortex current in the dual model by the for- 
mula, 



VxJ(t,r)= / dt' J d 3 r'Q(t-t',r-r')3 v (t',r'), (27) 
generalizing the classical static duality relation 

2tt/m,o J, = Vx J, (28) 

which is equivalent to Eq. (|2*3^l . since J = qp°v s . Note 
that the factor lix^ffi, accounts for an elementary flux 
quantum in the dual model corresponding to a unit 
charged Eq. I|28|) also follows from the continuum limit 
of the exact duality transformation on the lattice. 

The low-frequency limit of Q{ui) = Q{ui, p = 0) is given 

by 



= r oP °V 2 v s + qE, (22) 

where Tq is the bare kinetic coefficient, E is the electric 
field, and we have neglected the noise for simplicity. The 
superfluid velocity is related as usual to the phase of the 
order parameter, i.e., v s = V#. Because of the vortices, 
V x v s / 0. The vortex current is given by 



= I^Vxv, (23) 



Q(w) « ^g^fpl = 2np A . (29) 

Thus, when the fluctuations are taken into account, Eq. 
B28J1 will hold approximately in the low-frequency limit 
with pa,o replaced by pa- 

The electric-magnetic duality in the response functions 
implies a linear response 

J v (u,p)=a(v,p)B(u,p), (30) 
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where a is the dual AC conductivity. The above equation 
is a generalization of Eq. (|25[) . In order to check the 
validity of Eq. (|3()|> . we rewrite Eq. Ij26() in real space: 

J(t, r) = J dt' J d 3 r'a{t -t',r- r')E(t', r'), (31) 

where E = —dA/dt. By taking the curl of Eq. i|3T|) 
we obtain, after some trivial algebra, precisely Eq. (|30|l . 
with the dual conductivity related to the conductivity of 
the original model through the formula 



27T/M 



cr(w), 



(32) 



where we have assumed a low frequency regime at zero 
momentum. 

Since the DFT Hamiltonian (|20|l has the same form 
as the GL model in the MCF regime, we have that an 
anomalous dimension rjh — 1 is generated for the gauge 
field h. Therefore, the scaling behavior of the AC conduc- 
tivity in the MCF regime given by Eq. Ijl7|) also applies 
here and we obtain that <j(oj) ~ H 2 ~ z , where z denotes 
the dynamic critical exponent of the DFT. Thus, sim- 
ple dimensional analysis in Eq. (|32ll . using the scaling 
relation for d = 3, yields 



(33) 



Note that the dynamic exponent z drops out in the power 
counting. The above result would also follow by consid- 
ering a renormalized version of Eq. I|25|l . where ha,o is 
replaced by /i^. In the extreme type II limit ha ~ M 1 ^ 2 , 
which should have the same scaling as a ~ /i 2 ~ z , imply- 
ing once more that z = 3/2. 

The question now is how the exponent z is related to 
the exponent z of the dual model. We have seen that 
the static exponents coincide, for the obvious reason that 
both original and dual models lead to the same singular 
contribution for the free energy. The relation between z 
and z is much less obvious. It can be obtained as folows. 
From Eq. and the Joscphson relation we obtain that 
the right-hand side of Eq. (|32|l scales for d = 3 and p = 



er(w) ~ {-ilo) 1 ' 2 *. 



(34) 



Since a(u>) must scale like in the MCF regime, we can 
use Eq. (|12f) and write also 



(35) 



Comparison of Eqs. i|34[l and 13511 leads to the scaling 
relation 



Insertion of the exact result l|33(l in the above equation 
gives then the exact result 



(37) 



It is interesting to note that the scaling relation i|36|) 
still gives a reasonable result for z when the mean-field 
value z = 2 is used, leading to z = 8/5 = 1.6, which is 
not too discrepant with the exact value 

The above results imply the following duality relation 
between the conductivities of the original and dual mod- 
els, 



o~(uj)a{uj) ~ const, 



(38) 



as the critical point is approached. The above duality 
relation is analogous to the well known Dirac relation 
between electric and magnetic charges 



IV. DISCUSSION 

In this paper we have combined RG and duality ar- 
guments in order to solve a controversial issue, namely, 
on the value of the dynamic critical exponent z in ex- 
treme type II superconductors. Our analysis implies that 
z = 3/2 exactly, therefore confirming the result given by 
the Monte Carlo simulations of Ref. |(| and the exper- 
imental value of Ref. 0- An interesting consequence of 
our analysis is that the original model and its dual share 
the same dynamic exponent. Thus, z, a, and v are the 
same in both models, while r\ is not the sained 

It remains to discuss how our analysis fits in the classi- 
fication of dynamic models of Hohcnbcrg and Halperin. 4 
We have obtained a value of z identical to the one of 
model E, which is a model critical dynamics for the spin- 
wave modes of a superfluid. Since the spin-waves are 
decoupled from the vortex degrees of freedom^ and the 
duality analysis of the critical dynamics takes precisely 
the latter into account, model E should not be expected 
to govern the critical dynamics of an extreme type II 
superconductor at zero external field. The right model 
critical dynamics is actually given by a generalization of 
model A in which the dynamics of the vortex loops are 
taken into account. A simplified version of such a model 
was considered in the beginning of Section III-B, in Eqs. 
(P%)) and I23J- These equations describe respectively the 
dynamics of the supercurrent and the vortex current. An 
interesting physical consequence of such a dynamics is the 
electric-magnetic duality in the linear response of the su- 
percurrent and the vortex current. Such a point of view 
leads to a modified dynamic London equation, which is 
derived from the two following Maxwell equations gener- 
alized in such a way to include also the vortex current, 



(36) 



(39) 
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V x B = J 



<9E 
~dt' 



(40) 



The remaining Maxwell equations are not affected by the 
vortices, since V • J v =0. Due to this, the electric- 
magnetic duality holds only in the above two Maxwell 
equations, ft would only hold in all the four equations 
if open vortex lines were present, in which case magnetic 
monopoles would be attached to the vortex line ends. 
This is obviously not the case here since the U(l) group 
is not compact^ 

The electric-magnetic duality of the above equations 
corresponds to E — ► B, B — ► — E, and J — > — J„. We 
will derive the London equation using the approximation 
where 3 V has a very weak dependence on r, i.e., by using 
Eqs. lHU and (SI- Thus, by taking the curl of Eq. (U 
and taking into account Eqs. I|25|) . Q28)). and V ■ B = 0, 
we obtain 



d 2 B 



V 2 B - /4 B 



= 0. 



(41) 



Note that the above dynamic London equation contains a 



damping term. The general case is of course more com- 
plicate and Eqs. l|30*jl and ffify must be used instead. 
Remarkably, the duality approach allows us to derive an 
exact value for the dynamic exponent just by using scal- 
ing arguments. However, in order to derive the scaling 
functions for the AC conductivity a more explicit calcu- 
lation is necessary. This will be the subject of a forth- 
coming work^S 

We conclude by saying that although the dual dynam- 
ics studied here is different from the quantum casej«i 
the result l|38[) applies also there. We should be note, 
however, that in the quantum case the dynamics can be 
derived from the quantum Lagrangian, which is not case 
for a classical system near the critical temperature. 
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